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A. Determination of the diffusion length 


I. Introduction 


The diffusion length Z of thermal neutrons in a moderator is defined as the 
average distance thermal neutrons will diffuse in the moderator from the point 
where they become slowed down to thermal energies to the point where they 
get captured. The quantity ZL according to diffusion theory is related to the 
other fundamental constants of the moderator by the following equation: 
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L* (1) 
where A; = transport mean free path 
A = absorption mean free path. 

According to the open literature the diffusion length for water has earlier 
been experimentally determined by GAMERTSFELDER and GOLDHABER [1] and by 
BERTHELOT et al[2]. The authors (1) report a value L=3 cm+ 10% and 
(2) 2.77 + 0.04cm. As all natural neutron sources give fast neutrons, the method 
used by the authors cited above essentially has been the same, to bring about 
a differential neutron source. Since the determination (2) of LZ has been done 
with a Mn-detector and for a special case it seemed very interesting to make an 
investigation under more general conditions. 


il. The method 


The diffusion equation for thermal neutrons is simplified in the stationary 
case and when the slowing down density —0 to 


4 + 
Ag i = (2) 
where 9 = density of thermal neutrons. 
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The solution of this equation for a point source in an infinite medium is 


te 
ia 
o (r) = const. - (3) 


where r = distance from the source of thermal neutrons. 
To bring about a source of purely thermal neutrons a differential measure- 
ment is made. At first the spatial distribution in pure water of the thermal 


neutrons originating from the slowing down of fast neutrons from a Ra-«-Be- | 


source is studied. Then an absorber (Cd) of thermal neutrons which is spherically 
symmetric in relation to its source is introduced, and the spatial distribution 
of the thermal neutrons outside the absorber is studied. The difference between 
these two measurements represents the density of therma] neutrons captured 
by the absorber. This density will satisfy the equation (3). 

Throughout this investigation the term “thermal neutrons’ is used. As this 
term is used in a rather diffuse way in the literature it is worth while men- 
tioning that a better term covering the experimental conditions in this experi- 
ment would be C-neutrons, a term due to Fermi. In reality it is not possible 
to state that only thermal neutrons in the strict sense are measured in dif- 
ferential experiments of this type. The correct statement is that the energy of 
the measured neutrons is below the Cd-absorption limit. Below this limit, the 
energy spectrum of the neutrons is unknown and is in fact strongly disturbed 
by the absorbing cadmium medium, as will be shown by the measurements. 


If. Experimental arrangement 


Two neutron sources were used, 100 mC of Ra mixed with 500 mg Be and 
250 mC Ra mixed with 1250 mg Be, contained in cylinders of monelmetal, the 
external dimensions of which were 35 mm long, 14 mm wide, and 40 mm long, 


Fig. 1. The boron chamber. Total length 30 mm, diameter 7mm. Both electrodes of platinum. 
1. Collecting electrode. 2. High tension wire. 3. Cylinder. 4. Guard ring. 


Fig. 2. Mounting of the boron chamber at the end of the coppertube. 1. Polysterene tube. 


2. Coppertube. 3. Copperfoil. 4. Collecting wire. 5. High tension wire. 6. Polysterene plug. 
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to high tension 
to head amplifier 
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Fig. 3. Experimental arrangement. 1. Watertank with steel rails. 2. Carriages. 3. Steel 
scale. 4. Ra-%-Be-source (the spheres indicated by dots). 5. Boron chamber. 6. Pre-amplifier. 


17 mm wide, respectively. In order to bring forth a differential neutron source, 
two spheres were manufactured, each of 90 mm radius, one made of 0.4 mm 
Al-sheet, the other of 0.7 mm Cd-sheet. With a thin rod of polysterene, at one 
end furnished with a small cylindrical holder, the sources could be fixed with 
their centres in the centres of the spheres. The other end of the rod was rigidly 
connected to a carriage, movable along two iron rails over a U-shaped watertank 
of stainless steel, the dimensions of which were: 80 cm long, 40 cm wide, and 
60 cm deep. During the measurements this tank was filled with 140 liters of 
distilled water. As detector a small boron chamber was used, the construction 
of which is seen in Fig. 1. The chamber was filled to 1 atm. with BF;. It was 
fixed to one end of a thin coppertube which by a copperfoil was divided into 
two halves containing high tension and collecting wires, respectively. The other 
end of the tube was rigidly connected to another carriage which was also 
movable along the rails. Both carriages had a pionter showing directly on a 
steel scale the distance between the centre of the source and the centre of the 
boron chamber. The chamber as well as the copper tube were surrounded by 
a thin polysterene tube to prevent water from penetrating into the tube con- 
taining the wires (Fig. 2). Directly on the carriage holding the boron chamber 
a pre-amplifier was fixed from which cables led to the head amplifier (total 
amplification of the primary impulse ~ 10°) and a scale of 64 with discrimi- 
nator. The whole arrangement is shown in Fig. 3. With this apparatus two 
main series of measurements were made. First, the neutron density was measured 
with the Al-sphere filled with solutions of CdSO,4 in varying concentrations. 
Neutron source 100 mC Ra-«-Be. Secondly the density was studied outside the 
waterfilled Cd-sphere. Neutron source 250 mC Ra-«-Be. 
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IV. Results of the measurements 


A. Measurement with Al-sphere and CdSO, 


First the neutron density in pure water was measured. The same number off! 
pulses was counted either when the measurement was made with the waterfilled — 


Al-sphere or without it, which shows that the disturbance in neutron density | 
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Fig. 4. Neutron density curves measured in pure Fig. 5. The diffusion length as a func- 
water and outside the Al-sphere with Cd-solutions. tion of the ratio of concentration of 
The two upper curves give r-@(r) as a function of 7, cadmium to hydrogen in the sphere. 


caused by the thin sphere was negligeable. Then the sphere was filled with 10, 
20, 37, 50, 62, 100, 150, and 200 g CdSQ, dissolved in the 3 litres of water 
contained in the sphere. The result, represented by the group of curves, is 
shown in Fig. 4. If the difference between the density in pure water and the 
density outside the sphere, as measured point by point, is represented by o(r), 


it will be seen that 7-o0(7) = const. e ” according to theory shall form a line 
when plotted on semilogaritmic paper. In Fig. 4 this line is shown for the 
highest end lowest concentrations of Cd in the sphere. The values of Z for 
the various concentrations of Cd calculated by the method of least squares are 
given in table 1. 

L, thus, becomes a function of the Cd-concentration in the sphere as is 
graphically shown in Fig. 5. 

This dependence of L on the Cd-concentration must be intimately connected 
with the idealization made in the interpretation of the experimental results. 
Equation (3) predicts the neutron density distribution from a point source of 
thermal neutrons but in this experiment a sphere is substituted for the point 
source. As will be seen, there are at least two possibilities of making system- 
atical errors in these measurements. 
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Table 1 


Experimental values of Z for the different Cd-concentrations in the sphere: 
Coa and Cy atomic concentrations of Cd and H. 


Coa 
eg CdSO, ok 10* LT cm 
H 
10 146 3.20 + 0.07 
20 2.92 3.18+ 0.09 
37 5.41 3.16 + 0.09 
50 1.32 2.99 + 0.08 
| 62 9.07 | 2.84+0.07 
100 14.6 2.76 + 0.07 
| 150 22.0 2.80 + 0.06 
200 29.3 2.80 + 0.06 


a. In order that the measurements shall give the correct value of ZL, it is 
necessary that the distribution of sources of thermal neutrons outside the sphere 
at all points and at all Cd-concentrations is the same as the source distribution 
in pure water. Otherwise the density o(7) would be erroneous, since in taking 
the difference of the two densities with and without Cd-solution it is anticipated 
that nothing is altered in the neutron distribution except that a certain percent- 
age of thermal neutrons is absorbed in the sphere. 


b. Secondly, a great error may be caused by the large extension of the 
source. Neutrons slowed down to thermal energies in the neighbourhood of 
the sphere but outside of it have great chances to diffuse back into the sphere 
and get captured in Cd; in the case of pure water, on the other hand, these 
neutrons are able to diffuse back into the region where they can be counted. 
The sphere thus excerts a shadowing effect upon neutrons slowed down in its 
immediate neighbourhood. This effect must depend upon the solid angle occupied 
by the sphere as seen from the point of measurement. Thus, it contributes to 
bring about an erroneous reduction in neutron density in the vicinity of the 
sphere. It is to be expected that this error is considerable up to a distance of 
a diffusion length from the surface of the sphere and at high Cd-concentrations. 


B. Investigation of the source distribution at various conditions 


Since it is impossible to measure the source distribution for thermal neutrons 
directly, the only possibility to get an idea of this is to measure the distribu- 
tion of resonance neutrons with a detector of a material having a sharp re- 
sonance immediately above the thermal energy region. For this indium was 
chosen which has a pronounced resonance peak at 1.44 eV® and the absorption 
cross section of which is very high, ~ 10000 barns, so that a reasonably high 
activity is induced even at low neutron densities. With In-detectors (25 x 25 mm, 
48 mg/cm?) screened by 0.3 mm Cd-plates, the spatial distribution was measured 
outside the Al-sphere filled with pure water, 10 g CdSO, and 150 g CdSO, in 
solution. The result of these measurements is graphically shown in Fig. 6, 
As will be seen, the density of these neutrons at the higher Cd-concentration 
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Fig. 6. Density of In-resonance neutrons outside the Al-sphere. A. Pure water. B. 10 g 
CdSO, in the sphere. C. 150 g CdSOy, in the sphere. 


in the sphere differs markedly from the density in pure water. At r= 10 cm 
the density is lowered by some 10 %, which decrease does not vanish until 
r~17 cm. At the lower Cd-concentration the deviation is about 3 % at 
y=10cm and vanishes at r~14 cm. The reason for this decrease in the 
density of In-resonance neutrons might be that the absorption cross section of 
Cd is approximately 50 times the absorption cross section of H for neutron 


energies > 1.5 eV (=  64~ 91 for 2 ~ 9910-4 and ~0.03 for 2 15 10 
Ca Cy Cy Cy 


both relations giving eye 50), or else Cd might have less pronounced resonances 
OH 


for higher neutron energies. The perturbation in the source distribution de- 
creases with Cd-concentration. The result of this perturbation is that the ther- 
mal neutron density is more decreased near the sphere than at some distance 
from it, and thus g(r) is increased in the same manner which tends to give 
a too low value of L, especially at high Cd-concentrations. 


C. A study of the shadowing effect 


In order to get a value of LZ under different absorption conditions and, if 
possible, to show the second disturbance — the shadowing effect — more 
clearly, a series of measurements were made with the Cd-sphere. As in the 
first experiment (A), the density distribution of thermal neutrons was measured 
in pure water and outside the waterfilled Cd-sphere. The density curves are 

7 


shown in Fig. 7. The curve @(r):7 = const. e ", given in Fig. 8, has a peculiar 
appearance. At 7 ~ 13.5 cm the line is broken. For 10 <7 < 13.5 cm the line 
has a greater slope than for r > 13.5 cm. But this must be the case if the 
shadowing effect is considerable, because then the density in the neighbourhood 
of the sphere undergoes an extra decrease, which appears as an increase in 
o(r) and consequently a greater slope of the line. Furthermore, it is seen that 
this nick of the curve occurs at a distance from the sphere of ~ 45 mm, 
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Fig. 7. Neutron density in pure water and outside the Cd-sphere. 
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Fig. 8. The line giving ZL. The nick on the line occurs at approximately one diffusion 
length from the surface of the sphere. 


which approximately equals the diffusion length. In order to determine JL, 
only points for which r = 13.5 cm are used. The method of least squares gives 


ij—=3 29). 0.09 em 


a value which is in good agreement with the value 3.20 cm received at the 
lowest Cd-concentration in the first experiment (A). 

In order to further investigate the shadowing effect caused by the Cd-sphere 
another measurement was made. With empty spheres the neutron density was 
measured outside both the Al-sphere and the Cd-sphere. As is shown in Fig. 9, 
a very strong depression in neutron density is caused outside the Cd-sphere. 
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Fig. 9. Neutron density outside empty spheres. 


At r= 10 cm the decrease amounts to ~ 50 %, at r=13.5 cm it amounts to 
~10% and at r=18 cm it is ~1 %. In the case when the sphere is empty, 
however, the effect is much more pronounced than when the sphere contains a 
moderator, because all neutrons from the source in the former case are slowed 
down outside the sphere. Moreover, the absorption of the sphere is noticeable 
at larger distances from it, owing to the much higher density of thermal neu- 
trons. In any case it will be seen that the source of error caused by the 
shadowing effect must be taken into consideration, when the goal is to deter- 
mine L by this method. The best way to eliminate the two errors now men- 
tioned, is to make the measurements at larger distances from the sphere. But 
at 7 > 20 cm the difference @(r) becomes very small and the statistical errors 
make a measurement impossible. Therefore, it seems that the best value of L 
is obtained if points are used for which 13.5 =7r< 19.5 em. 


V. Discussion of the results 


Now that the two sources of error have been investigated, it is possible to 
discuss the curve in figure 5. At higher Cd-concentrations in the sphere prac- 
tically all neutrons slowed down outside of it that happen to diffuse back into 
the sphere will be absorbed, and the perturbation in thermal sources is rather 
constant. In that region the Z-value will be constant. At a certain concentra- 
: ee : : ; 
tion ratio, a orn 9.10~-4, the shadowing effect will begin to decrease because 

H 
some thermal neutrons can diffuse back out of the sphere and get counted, 
which means that o(r) decreases comparatively much in the neighbourhood of 
, Ce : 
the sphere and the value of Z increases. When ae ~5-10~4, the absorption 
H 
probability in Cd then being only 5 times the absorption probability in H, 
the neutrons are able to diffuse to and from through the wall of the sphere 
and have a greater chance of being counted. Thus, the value of Z will increase 


only slowly with the Cd-Concentration for ee <5-10-*. Even at low con- 
H 


54 


ARKIV FOR FYSIK. Bd 2 nr 7 


centrations of CdSQ,4, there remains a disturbance in the source distribution. 
The measurement with the Cd-sphere, where this disturbance is inappreciable 
for r>13 cm, gives a value of ZL which must be very probable on the as- 
sumption that the neutron density is measured at such large distances from 
the sphere that the shadowing effect is inappreciable or eliminated. The mea- 
=), which must be an- 
Cy 


other very probable value, since all perturbations in neutron density vanish 
: C 

with Cca. The curve L=f (ce 

H 

giving L = 3.22 cm. 


surement with the Al-sphere gives the limit lim L( 
Coq>0 


is for that reason extrapolated to Coq = 0 


VI. Discussion of the other experimental errors 


With the value of Z measured as outlined above one may assume that the 
two greatest systematic errors are eliminated or very small. 

A third systematic error might be caused by the finite size of the boron 
chamber. This would tend to give too high neutron densities. On the other 
hand, a small fraction of the neutrons must get absorbed in the glass walls 
of the chamber. Moreover, a calculation shows that the absorption of thermal 
neutrons by the boron in the chamber is equivalent to the absorption in the 
corresponding water volume (Ny ox ~ Ngosp). For these reasons the error caused 
by the finite size of the detector must be very small. No correction is applied 
for it. 

The counting of the boron chamber, due to neutrons of energies greater 
than the cadmium absorption limit, was <1 % of the total counting. More- 
over, since the distribution of non-thermal neutrons is little disturbed by in- 
troducing the absorbers, the error introduced by counting non-thermal neutrons 
must be negligible. 

The statistical errors were kept very small, since least 10000 pulses were 
always counted, which corresponds to a statistical error of 0.5—1 %. The 
background count was 0.04 impulses/min. The statistical uncertainty in 0 (r) 
thus amounts to 2 % at r=10.0 cm and to 5 % at r=18.5 cm. During the 
whole series of measurement the reproducibility of counting was checked and 
all measurements were made several times. The measured number of pulses 
always fell within the statistical limits of error, which guarantees the reliability 
of the electronic equipment. 

The errors in 7 amounts to + 0.4 mm, the relative error in r thus amounts 
to 0.4 % at most. 

The theory of this experiment requires a slowing down medium of infinite 
size, a condition which is not satisfied. When, however, only thermal neutrons 
are studied and the detector never approaches the walls of the tank closer 
than six diffusion lengths, this disturbance in boundary conditions is neglible. 

The total uncertainty in L is estimated to be 4 %. 

The measurements described above then give as a weighted mean 
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which means, correctly speaking, that the diffusion length of C-neutrons under 
the experimental conditions specified above is 3.25 cm, this value being cor- 
rected for the disturbances in the energy distribution of the C-neutrons that 
arise during measurements. This value thus agrees with the value reported 
by GAMERTSFELDER and GoLDHABER[1] within limits of error. The value 
2.77 + 0.04 cm reported by BeErrHELor et al[2] agrees with the value 
2.80 + 0.06 found in the measurements with the Al-sphere at higher Cd- 


concentrations. This value, however, is not corrected for the two great disturb- — 


ances in the neutron distribution as shown in these measurements. 
Furthermore the diffusion theory approximation suffers from a theoretical 


ae : . Oa : : 
error, which is a function of the ratio — (oq = absorption cross-section, ¢ = total 
o 


cross-section). [8] This error might be estimated to ~ 1 % in the mathematical 
description of the density of neutrons in water. A better approximation in 
place of eq. 1, for instance, is 


4 hy 
lhe =A. u (1+ 4), 


this expression resulting from the general Bolzmann equation. [8] This theoretical 
error in ordinary diffusion theory of course increases the uncertainty in the 
value of L given above. 

B. Determination of the second spatial moment 


I. Introduction and theory 


By age-diffusion theory [4,5] the second spatial moment <r?> for thermal 
neutrons is functionally related to the diffusion length by the following relation 


rin? 


6 = Mi, = Tin + ive (4) 
where M;, = migration area 
Tin = ‘age’ of thermal neutrons 


L = diffusion length of thermal neutrons. 


In spite of the fact that age-diffusion theory is a poor approximation in 
problems of slowing down in a hydrogenous moderator like water — especially 
when the neutrons have so high initial energies as Ra-«-Be-neutrons — it 
seemed interesting to determine (rj). In fact this measurement with a small 
boron chamber might give an interesting comparation with results received in 
measurements with foils by several authors. [6, 7]. 

Experimentally the density @(r) of thermal neutrons is measured as a func- 
tion of the distance r from the source. The total number of neutrons counted 
by the boron chamber is then 


N=42 ee rdr 
0 
where J is the number of counts/min at r. 
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The mean square of the migration distance is then found from 


[Irtdr 

Tin) = 6 Min = = a 
{Irdr 
0 


II. Experimental equipment and results 


The whole apparatus used to determine <7?) is the same as described above 
in the determination of LZ. The density of thermal neutrons was measured 
from r= 4 cm to r= 38 cm. For 4<r < 20 cm the 100 mC Ra-a-Be-source 
was used. For 12 =r <38 cm the 250 mC neutron source was used in order 
to get reasonably high counting rates. In the region 12 < 7 < 20 cm the two 
sources were calibrated against each other. The result of this calibration is 
given in the table. 


Table 2 


Calibration of the 250 mC source against the 100 mC source. 


rote I counts/min I counts/min Li00 
100 mC 250 mC Teen 

12.0 947 2S 0.436 

13.0 758 1720 0.440 

14.0 611 1385 0.441 

15.0 | 490 1123 0.436 

16.0 398 911 0.438 

| rie) 321 734 0.437 
18.0 269 614 0.438 
19.0 216 490 0.441 

20.0 | 179 415 0.432 

| Mean value 0.438 


The readings for r > 20 cm with the 250 mC source were adjusted to the 
proper values corresponding the 100 mC source by multiplication with the 
factor 0.438. The electronic apparatuses run continously during the course of 
the experiment, and no drift could be detected within the statistical accuracy 
of the counts. At each point up to r= 25 cm 10000 to 20000 pulses were 
counted. The least number of pulses counted for 7 > 25 cm was actually 
7000. The statistical accuracy thus amounts to ~1 %. In Fig. 10 the function 
Ir? =f (r) is shown on semilogarithmic scale. This curve shows an exponential 
decrease when 7 > 20 em with the relaxation length 9.93 cm. Assuming an 
exponential decrease without limit with this relaxation length, points on the 
curve Ir? = f(r) were calculated for r > 38 cm and J r* was formed. However, 
Ir* = y(r) shows no clear exponential decrease until r > 86 cm (Fig. 11). The 
space integral (Ir?dr) was evaluated by formal integration for r > 38 em, 
assuming exponential decrease without limit, and by Simpsons rule over the 
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Fig. 10. Semilog plot of I7? =f (r), the curve showing exponential decrease from r = 20 cm. 
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Fig. 11. Semilog plot of Ir*= g(r). The dotted part of the curve is calculated from the 
exponential decrease of Ir? = f(r). The curve I rt = (r) shows no safe exponential decrease 
until 7 = 86 ecm. 


remainder. (Ir*dr) was evaluated by Simpsons rule in the region 0 < 7 < 86 cm 
and by formal integration over the remainder. The result is: 


oth 973.4 - 10° 
Ne 0 = Or ae ; 2 
Cin) = 7 © Prong pr poe ei 
bap LF 
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The corresponding migration area is 


2 
MP ce Se = 61.0 cm?. 


The above value of Z inserted in equation 4 gives 
Tin = (61.0—10.8) cm? = 50.2 cm?. 


This value of <7j,> is some 10 % higher than the value reported by Rusu [6] 
and by Munn and Pontscorvo [7]. As is easily seen there are to be expected 
relatively large discrepances between this measurement with boron chamber 
and the measurements [6,7] with In- and Dy-detectors respectively. With Dy- 
detectors, for instance, not only the thermal neutron distribution is measured, 
but the distribution of resonance neutrons is superposed, because the cross 
section of Dy shows disturbing resonances at a few eV°. This tends to give 
lower values of <7rj,)>. The same reasoning holds for In-detectors. Boron, on 


: Reet | : ; 
the other hand, is a strict > absorber practically over the whole energy region 


of the Ra-«-Be-neutrons. As measurements show, the counting rate of the 
Cd-shielded BFjchamber is <1 % of the unshielded chamber, so that the 
error introduced by counting non-thermal neutrons is very small in the case 
of a boron-chamber. 
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